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• Abstract 

^ . Motivated by the data from Super-Kamiokande and elsewhere indicating oscihations 

H ' of atmospheric and solar neutrinos, we study charged-lepton-flavour violation, in 

■ - - ' particular the radiative decays fi ^ and r — >■ /i7, but also commenting on /i ^ 3e 

and r — > 3/i/e decays, as well as — > e conversion on nuclei. We first show how the 
renormalization group may be used to calculate flavour- violating soft supersymmetry- 
breaking masses for charged sleptons and sneutrinos in models with universal input 
parameters. Subsequently, we classify possible patterns of lepton-flavour violation in 
the context of phenomenological neutrino mass textures that accommodate the Super- 
Kamiokande data, giving examples based on Abelian flavour symmetries. Then we 
calculate in these examples rates for /i — > 67 and r — > which may be close to 
the present experimental upper limits, and show how they may distinguish between 
the different generic mixing patterns. The rates are promisingly large when the 
soft supersymmetry-breaking mass parameters are chosen to be consistent with the 
cosmological relic-density constraints. In addition, we discuss fj, ^ e conversion on 
Titanium, which may also be accessible to future experiments. 
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1. Introduction 



There has been increasing interest in massive neutrinos during the past year, triggered 
principally by the Super-Kamiokande data on the z/^/z/g ratio in the atmosphere. The 
latter was found to be significantly smaller than the Standard Model expectations, with a 
characteristic azimuthal-angle dependence indicating the presence of neutrino oscillations. 
The data analysis favours —>■ Vj. oscillations, with parameters in the ranges 

- (10-2 to 10-3) eV^ (1) 

sm^2e^r > 0.8 (2) 

Dominance by oscillations is disfavoured by the Super-Kamiokande data on 

electron-like events as well as by the data from the Chooz reactor experiment 0. 
Oscillations involving a sterile neutrino are disfavoured, but not yet excluded, by a detailed 
study of the azimuthal-angle dependence of muon-like events [|l| and by measurements of 
7r° production. Moreover, in most theoretical models sterile neutrinos tend to be heavy. 
Therefore, we consider — > z/,- as the 'established' hypothesis for the atmospheric neutrino 
data. 

In addition, the long-standing deficit of solar measured on Earth may also be ex- 
plained via neutrino oscillations, either in vacuo or enhanced in matter by the Mikheyev- 
Smirnov-Wolfenstein (MSW) mechanism [Q. The first option would require 6ml^^,^ ~ 
(0.5 - 1.1) X 10-1° eV^ where q; is /i or r. MSW oscillations [^, on the other hand, require 
Sml^^^ ~ (0.3 - 20) X 10"^ eV^ with either large sin^ 2^ea ~ 1 or small sm'^2e^a ~ 10"^. 
The presence of either z/g or v^, —>■ oscillations at a high level is, therefore, an open 
question. 

Both the solar and atmospheric neutrino data can be accommodated in a natural way 
in schemes with three light neutrinos with at least one large mixing angle and hierarchical 
masses, of the order of the required mass differences: ~ (10-^ to lO-^'^) eV and 
777,2 ~ (10"^ to 10"^) eV ^ m^. On the other hand, if neutrinos were also to provide 
significant hot dark matter, three almost-degenerate neutrinos with masses of ~ 1 eV 
would be needed. 

Neutrino oscillations involve violations of the individual lepton numbers -Z^e,/i,T; raising 
the prospect that there might also exist observable processes that violate charged-lepton 
number conservation 0, such as /i ^ 67, 3e and r — > fi'y,3fi/e, and /i — > e conversion 
on heavy nuclei |], 0, ^ |l9. We recall briefiy the present experimental upper limits 
on the most interesting of these decays for our subsequent discussion: 

5/?(/i 67) < 1.2 X 10-11 : [11] (3) 
BR{i2+ e+e+e-) < 1.0 x lO-^^ : ^ (4) 
R{l^~Tt e-Ti) < 6.1 X lO-^^ : O (5) 



BR{t fi-f) < 1.1 X : (6) 



-2- 



Projects are currently underway to improve several of these upper limits significantly: 

BR{i2 67) 10"^^ : [15] (7) 

R{lJ,-Ti e'Ti) ^ few x 10"^^ : |T6l (8) 

R{lJ,-Al e-Al) 1 X 10"^^ : 

BR{t ^ /i7) ^ 1 X 10~^ ? : 



(9) 

(10) 



and there are active discussions of intense /x sources that might enable the upper limits on 
/i — i> e transitions to be improved by several further orders of magnitude \V5\. 

We evaluate the possibility of charged-lepton-fiavour violation using the most natural 
mechanism for obtaining neutrino masses in the sub-eV range, namely the see-saw mecha- 
nism 1 2^, which involves Dirac neutrino masses of the same order as the charged- lepton 
and quark masses, and heavy Majorana masses Mj,^, leading to light effective neutrino mass 
matrices: 

m,ff = m^-{M,^)-^-mf. (11) 

Neutrino-flavour mixing may then occur through either the Dirac and/or the Majorana 
mass matrices, which may also feed flavour violation through to the charged leptons. 

The specific mechanism explored in this paper is renormalization of the sneutrino and 
slepton masses in a supersymmetric theory via the neutrino Dirac couplings P|. It is 
well known that the prototypical charged-lepton flavour-violating process — > 67 provides 
one of the most stringent upper limits on flavour violation in the Minimal Supersymmetric 
extension of the Standard Model (MSSM). If the soft supersymmetry-breaking sneutrino 
and slepton masses were non-universal before renormalization, it would be very difficult to 
understand why this decay was not seen long ago. Universal scalar masses arise naturally 
in no-scale supergravity models pl[, the framework favoured here, as well as in gauge- 



mediated models ||2^ . In the universal supergravity case, the soft supersymmetry-breaking 
sneutrino and slepton masses are subject to calculable and non-trivial renormalization via 
Dirac neutrino couplings at scales between Mqut ~ 10^^ GeV and M^^ ~ lO^'^ GeV. 

The predictions of this class of universal supergravity models are quite characteristic. In 
non-supersymmetric models with massive neutrinos, the amplitudes for the charged-lepton- 
fiavour violation are proportional to inverse powers of the right-handed neutrino mass scale 
0] • Since the latter is much higher than the electroweak scale, the rates for rare decays 
such as /i — > 67 are extremely suppressed @]. On the other hand, in supersymmetric models 
these processes are only suppressed by inverse powers of the supersymmetry breaking scale, 
which is at most 1 TeV Among such models, those with non- universal input scalar 
masses at the GUT scale generally predict excessive rates for rare charged-lepton-fiavour 



violation, whereas they are very suppressed in no-scale pTf and gauge-mediated models [|2 



The class of supergravity models with universal scalar masses that we consider here toe the 
fine line between excessive and unobservable charged-lepton-fiavour violation, as we discuss 
in more detail below. 

Here we re-analyze the prospects for charged-lepton-fiavour violation in this theoretical 
framework, using as a guide the indications from Super-Kamiokande and elsewhere on 
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neutrino masses and mixing. Within the general see-saw scenario, solutions with various 
differences in the neutrino mass-matrix structure have been proposed, including models 
with maximal ||23|, , close-to- maximal |2^, ^ and bi- maximal neutrino mixing. In 
the following, we categorize models according to whether the off-diagonal elements in their 
Dirac and Majorana couplings 'match' in such a way that their mixing is almost two- 
generational (and hence predominantly in the yU — r sector), and more general models in 
which they 'mismatch', and the mixing is essentially three-generational (and substantial 
also in the fi — e and r — e sectors). We provide examples of the 'matched' and 'mismatched' 
categories in schemes with Abelian pS, ESl |2H] flavour symmetries, and comment on the 



possibilities with non- Abelian p9| flavour symmetries, as well as in a string-derived flipped 
SU{5) model whose characteristic property is the appearance of large off-diagonal /i — e 
couplings in the Dirac neutrino-mass matrix [^. As we exemplify with calculations in the 
Abelian models, the rates for the radiative decays /i — 67 and r ^ /i7 may offer good 
prospects for testing different textures. These decays may well take place at observable 
rates, and different neutrino-mass models correlate their decay rates in characteristically 
different ways, enabling —>■ ej and t ^ fij to serve as useful diagnostic tools for neutrino- 
mass models. Finally, we calculate in some models the rate for fi ^ e conversion on 
Titanium, which may also be accessible to future experiments, and comment briefly on 
/i ^ 3e decay. 



2. General Aspects of Charged-Lepton-Flavour Violation in Su- 
per symmetric Models with Universal Breaking 

We first display in Fig. 1 the one-loop diagrams that give rise to /i — >■ 67, noting that the 
T /i7-decay is generated by an analogous set of graphs. We later extend the discussion 
to include /i — 3e decay, r 3/i/e and /i ^ e conversion. The matrix element of the 
electromagnetic-current operator between two distinct lepton states k and Ij is given in 
general by 

Tx = {h\{p - q)\JxUp)) 

= Ui{p- q){m,iaxf3q^ {^mPl + ^f/^^i?) + 

(g%A - qxl ■ q) {A^ePl + A^Pr) }u,{p) (12) 

where q is the photon momentum. The coefficients Am and Ae receive contributions 
from both neutralino (n) /charged-slepton (Fig. 1(a)) and chargino (c)/sneutrino (Fig. 1(b) 
exchanges: 

aL^R aL^R I aL^R aL^R aL^R I aL^R /-I q\ 

— ^M{n) + ^M(c)' — ^E{n) + ^E(c) y^'^l 

The amplitude of the process is then proportional to T\t^, where is the photon-polarization 
vector. 

We recall that the easiest way to determine the loop momentum-integral contributions 
to the coefficients Am,e is to identify, in the corresponding diagram, terms proportional to 
{p ■ e) and {q ■ e). The coefficient of the former is proportional to the momentum-integral 



-4- 




1 1 



LI ' ^ e Li ' 
i_ I 1 JZ I— 




5^(„) ^ic) 
(a) (b) 



Figure 1: Generic Feynman diagrams for fi ^ cj decay: I represents a charged slepton (a) 
or sneutrino (h), and x*-"-* and x'^^'> represent neutralinos and charginos respectively. 



contribution to the axp term in (|T2]), while the coefficient of the latter is proportional to 
the difference between the momentum-integral contributions to the axp and (g^7A — qxl ■ q) 
terms. Defining x = M'^/m'^, where M is the chargino (neutralino) mass and m the 
sneutrino (charged slepton) mass, the following functions appear in the Am terms [ffn, H]: 



^M(n) '■ .\4 (1 — 6x + 3x^ + 2x^ — logx) and 



jT^il -x^ + 2x logx)^ 
^M{c) '■ (2 + 3x — 6x^ + x^ + 6x log x) and 

jr^(-3 + Ax - x^ - 2 log x)^ 

where mi. is the mass of the Ij lepton, while for the Ae terms we have: 

AE{n) ■■ (2 - 9x + 18a;2 - 18x^ + Qx^ log x) 

Ae{c) ■■ (13^(16 - 45x + 36x^ - 7x^ + 6(2 - 3x) logo;) 



(14) 



(15) 



Note in this case the lack of terms proportional to the gaugino mass M. The branching 
ratio (BR) of the decay Ij ^ /j + 7 is then given by: 

BR{1, ^ la) = ^ {{A^M? + (16) 

We see, therefore, that the branching ratios for radiative lepton decays involve the masses 
of several supersymmetric particles at low energies. 

As stated in the Introduction, in this work we assume universal scalar masses and tri- 
linear terms A at the GUT scale. However, the physical values of these masses to be used 
in (|l^,^,0) have to be obtained by integrating the renormalization-group equations of the 
MSSM supplemented with right-handed neutrinos, found, for example, in |^). The Dirac 
neutrino and charged lepton Yukawa couplings cannot, in general, be diagonalized simul- 
taneously. Since both these sets of lepton Yukawa couplings appear in the renormalisation- 
group equations, the lepton Yukawa matrices and the slepton mass matrices can not be 
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simultaneously diagonalized at low energies either. Indeed, in the basis where me is diago- 
nal, the slepton-mass matrix acquires non-diagonal contributions from renormalization at 
scales below Mqut, of the form p: 



loTT^ Mat ' 

where A^) is the Dirac neutrino Yukawa coupling, M^v is the intermediate scale where the 
effective neutrino-mass operator is formed, and a is related to the trilinear mass parameter: 
Ai = am3/2, where my 2 is the common value of the scalar masses at the GUT scale. As a 
result, the diagrams of Fig. 1 lead to radiative decays of charged leptons Q. 

We obtain the physical charged-slepton masses by numerical diagonalization of the 
following 6x6 matrix: 

where all tha entries are 3x3 matrices in flavour space. Using the superfield basis where 
Xi is diagonal, it is convenient for later use to write the 3x3 entries of ([T8|) in the form: 



^II = {mjf + Sml^+m^^- Ml{-- sin^9w)cos2p (19) 



m|j^ = (ml-^y + mj- M^sin^ewC0s2P (20) 

= {Al + 6Ae - iJ.tan(3)mg (21) 

^LR = ^RL (22) 

where tan f3 is the standard ratio of the two MSSM Higgs vevs, (^«|)^, {ml-^^ and A^ denote 
the diagonal contributions to the corresponding matrices, obtained by numerical integration 
of the renormalization-group equations, and Sml^ and 6A1 denote the off-diagonal terms 
that appear because Xd and Xi may not be diagonalized simultaneously ([T7|) . 

The full mass matrix for left- and right-handed sneutrinos has a 12 x 12 structure, given 
in terms of 3 x 3 Dirac, Majorana and sneutrino mass matrices. The effective 3x3 mass- 
squared matrix for the left-handed sneutrinos has the same form as the mj^j^ part ( P^D of 
the 6x6 charged-slepton matrix (|18]), with the difference that now the Dirac masses are 
absent. One might have expected that - as in the case of charged sleptons - the Dirac terms 
would induce considerable mixing effects. However, we show here that this is not the case 
in the sneutrino-mass matrix. Due to the vastly different scales involved in the full 12 x 12 
sneutrino mass matrix, a complete analysis is not straightforward. Instead, we construct 
an effective 6x6 matrix for the light sector, by applying matrix perturbation theory, in 



"'^We note that a complete renormalization-group analysis would also involve the quark and squark 
sector, and require a treatment of supersymmetric thresholds. However, the inclusion of such detailed 
effects would not affect the main conclusions of our analysis, so we may neglect them for our purposes. 



-6- 



analogy with to the see-saw mechanism ||20|. To second order, the result is: 



i^Deff 



( m\ + 0{mjM-^) ((2A, + A;v) -2/icot/?)- \ 
((2A, + A^) -2^cot/5)- ml + 0{m\M-^) 



(23) 



The first- and second-order terms in (|2^) are obtained assuming that all the parameters 
as real and the matrices A^^tv are proportional to the identity. Notice that the second- 
order terms along the diagonal can be neglected, whereas the first-order off-diagonal terms 
must be retained, since they lead to complete mixing of the pair-wise degenerate states. 
However, this does not affect the branching ratios for the flavour- violating radiative decays. 
Therefore, we simply use 



ml = (mfY + 5ml^ + -Mlcos2p (24) 



in our subsequent calculations. 

Next we consider the rare decay /i 3e and fi e conversion on nuclei, starting with 
the fi 3e reaction. This reaction is interesting on its own since it has a structure much 
richer than that of the fi —>■ e'j decay. Thus, it is possible that n —>■ 3e can take place in 
cases where fi —>■ ej is forbidden. This and the related r — >■ 3/x/e decay receive contributions 
from three types of Feynman diagrams. The first are photon 'penguin' diagrams related 
to the diagrams for /i — 67 and r /i/e7 discussed above, where now the photon is 
virtual and decays into an e+e~ (or fj.'^fi~) pair. A second class of diagrams is obtained by 
replacing the photon line with a Z boson. Finally, there are also box diagrams. In addition, 
all the above types of diagrams are accompanied by their supersymmetric analogues. We 
evaluate all the relevant diagrams exactly in our subsequent numerical analysis. However, 
we note that the dominant diagrams in the models of interest to us here are generally the 
'penguin' diagrams with an intermediate off-shell photon, which contribute via the Am and 
Ae terms presented in (|l^|15D. Compared to /i — >• 67 decay, the branching ratio is, 

-— ^ ^ 6 X 10"^ 25 

r(^+^e+7) ^ ^ 

This does not necessarily mean that /i ^ 3e decay is uninteresting to experiment, because 
the experimental detection and background problems are very different for the two decays. 
However, we do not present detailed numerical results for /i"*" — >■ e^e^e~ decay, because the 
factor (|25| ) is essentially universal. There is a similar small ratio for r — > 3/i/e relative to 
r —>■ fJ'/e'j, which seems to preclude its observation even at the LHC 



The /i — i> e conversion is a coherent process in a muonic atom originally studied in [p2 



33| . Even though this reaction proceeds via the same classes of diagrams as those discussed 
above in connection with /i+ e~^e~^e~, this reaction is rather different, since it involves 
hadronic currents. Effects of nuclear nature, such as the size of the nucleus in particular 
when heavy atoms are involved, play important role. Detailed calculations including all 
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contributions from penguin and the box diagrams will be given in section 5. Here, in 
order to obtain a first rough estimate of the ratio of the /i — e conversion to the /i — 67 
reaction, we restrict to the photonic contribution which dominates over a large portion of 
the parameter space. The fi e conversion rate relative to conventional muon capture, is 
given by 

r(/i ^ Uf,) Gf mf^ CZf{A, Z) 

where A^/^ = A^^^^+A^^ . It is important to note that the combination of matrix elements 
A^^^ in p6|) is different from that in appearing /i — >• e decay (p!6D. The function f{A, Z) 
is given in IQ, and has the following approximate value for elements with A ^ 2Z 



f(A, Z) = 1.0 - O.OSA + 0,25 (A - 1) + 3.24 (A - 1 - ±| A _ i| , « 0.I6 (27) 



Further, the parameter C in (pBf ) is ||3^ 



C= |F^ + 3F] + F|-2FaFp| -5.9 (28) 

and the nuclear matrix element \ME\ is 

\ME\ = I < (A^)/|J(0)|(A,Z),, > I ^ ZF,(g2), (29) 

where, in the case of the photonic diagram, Fc{q^) stands for the proton or neutron form 
factor. Calculations using the nuclear shell model for ^^22Ti have yielded Fc{q^ = —fn'^^) ^ 



0.543 for the case of proton and 0.528 for the neutron [S^]. Comparing now with 
B{fi~^ —>■ e"^7), one obtains a rough estimate of the expected range of the fi e conversion 
ratio: 

371 mf^ Cf{A,Z) 
^ 5.6 X 10-^BR{fx 67) (30) 

which shows a relative suppression about two orders of magnitude. Nevertheless, /i ^ e 
conversion is also interesting. Current experimental bounds give R{fi'^Ti — > e^Ti) < 
6.1 X 10~^^, whilst ongoing experiments will reach ~ 2 x 10^^"^. However, with an intense 
proton (and muon) source, such as that projected for a neutrino factory or a muon collider, 
experiments sensitive to rates as low as 10~^^ may be feasible |]19|. Moreover, as is apparent 



from (^6[) above, and the appearance of box diagrams, the structure of the matrix element 



is different from that for n ^ ej (p!6|), so the ratio ( pO| ) is not universal, unlike the /i — >• 3e 
case. Therefore, although we emphasize fi —>■ ej in what follows, we shall also present later 
some results for fi —>■ e conversion on Titanium. 
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3. Neutrino Mass Textures in the Light of Super-Kamiokande 



Having discussed the theoretical framework for calculating charged-lepton-flavour vio- 
lation in a supersymmetric model with universal input soft supersymmetry-breaking pa- 
rameters at the GUT scale, we now discuss possible extreme patterns of neutrino masses, 
mixings and Yukawa couplings that might figure in the renormalization-group equations. 

The most basic piece of information from the recent atmospheric-neutrino data is the 
existence of at least one large mixing angle in the lepton sector, that associated with the 
— r flavour mixing. Lepton mixing may in general arise either from the charged-lepton 
sector, or the neutrino sector, or both. In analogy to the quark mixing matrix Vckm, the 



leptonic mixing matrix Vmns is defined as |]37 | 



Vmns = V,Vj (31) 

where Vi transforms the left-handed charged leptons to a diagonal mass basis, whereas K 
diagonalizes the light-neutrino mass matrix frieff. In the see-saw framework which provides 



a natural mechanism for generating very light neutrinos |20], the latter is given by 

m,ff = m^-{M,^)-'-mf. (32) 

where and M^^ stand for the Dirac and the heavy Majorana neutrino mass matrices 
respectively. 

How may one characterize the structures of Dirac and heavy Majorana matrices that 



generate viable neutrino textures consistent with the Super-Kamiokande data? In ||2J], we 
proposed a classification according to the criteria of 'matched' and 'mismatched mixing', 
as defined below. 

• Matched mixing: This occurs when there is only one large neutrino mixing angle, 
namely that in the (2-3) sector of the light neutrino mass matrix rrieff, as suggested by 
the atmospheric neutrino data, and there is no large mixing in other sectors of either 
the light Majorana or the Dirac neutrino mass matrix. In this case, the problem reduces 
approximately to a 2 x 2 mixing problem. Since there is no way to render three degenerate 
neutrinos consistent with the bounds from neutrino less double beta decay without also 
large (1-2) neutrino mixing, 'matched mixing' requires hierarchical neutrino masses. 

In this case, it has been shown that (in the absence of zero-determinant solutions, i.e., 
solutions where strong cancellations in the (2,3) sub-determinant in the neutrino sector 
cause one of the eigenvalues to be relatively small, which we will discuss in an example be- 
low) the lepton mixing originates entirely from the Dirac mass matrices, while the structure 
of the heavy Majorana mass matrix M^^ does not affect the low-energy lepton mixing [Q. 
This result remains valid even when renormalisation group effects are taken into account 
38| . Then, in the basis where the charged lepton flavours are diagonal, the Dirac neutrino 



mass matrix takes a very simple form, given by 




oc 



/ \ 

xy X 
\0 y I J 



(33) 
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for symmetric and asymmetric textures, respectively 0. 

It is evident that, in such a scenario, the (1-2) mixing and (1-3) mixing are both zero, 
in first approximation, so the large (2-3) mixing is not communicated at all to the (1-2) 
sector, for which even the sub-dominant contributions are very small. Passing to the (2-3) 
mixing, we see that it can be either maximal or non-maximal, depending on the value 
of X, which can be as large as unity. It is clear that such a 'matched- mixing' scenario is 
consistent only with the small-mixing-angle (SMA) MSW solution for solar neutrinos, since 
the large-mixing-angle (LMA) MSW solution and the vacuum-oscillation (VO) scenarios 
both require large (1-2) and/or (1-3) mixing. 

If (1-2) and (1-3) mixing are both small, as in a generic 'matched-mixing' model, the 
fi —>■ ej and r — * 67 rates should be relatively 'small', whereas the t —>■ fi'-f rate may 
be relatively 'large'. These general expectations are borne out in the model calculations 
presented later. 

• Mismatched mixing: Entirely different structures arise when (i) there is more than one 
large mixing angle in rrieff, and/or (ii) there is a large Dirac mixing angle that involves dif- 
ferent generations from those of the light Majorana matrix. A mild example of mismatched 
mixing occurs, for example, when the atmospheric problem is solved by z/^ —>■ Vr oscilla- 
tions, whilst the Dirac mass matrix is related to the quark mass matrix, with Cabibbo-size 
mixing between the first and second generations. 

In generic mismatched mixing models, there are relatively large violations of charged- 
lepton flavour in all the (1-2), (2-3) and (3-1) channels. Thus, such a 'mismatched-mixing' 
scenario is a priori compatible with either the LMA or VO solutions of the solar-neutrino 
problem. Moreover, /i — > 67 and r ^ 67 generically have larger rates than they would 
have in matched-mixing models. However, the structure of the Majorana matrix becomes 
more complicated in the case of mismatched mixing. In particular, it is possible that the 
Dirac mass matrix is almost diagonal, with a large hierarchy of Dirac couplings, so that, in 
particular, Ai <^ A2 where the Aj are the eigenvalues of the neutrino Dirac coupling matrix. 
In this case, the light entry of the heavy Majorana mass matrix again effectively decouples 
from the heavier ones However, this is no longer true if the (1-2) mixing angle in the 



Dirac mass matrix increases. 

One particularly interesting example of such large mixing, which has been extensively 
studied in the literature, is that of bi-maximal mixing for which one can easily obtain 



viable solutions with degenerate neutrinos. As an example, we quote the texture 

(34) 



me// oc 



1 1 _ 1 

V2 2 2 
\ 2 2 



^Here we have correlated the (2-2) with the (2-3) and (3-2) elements, assuming that they arise from a 
flavour symmetry of the type discussed later. 
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For this texture, in the neutrino mixing parametrization 



I 



C2C3 



C2S3 



S2e 



-C1S3 - SiS2C3e*'' +C1C3 - SiS2S3e*'' 

\ +S1S3 - cxs^c^e'^ -S1C3 - CiS2S3e*'' 




(35) 



where the diagonal matrix ^g/T diag {rriie^'^, 77126^'^' ,771^), and and 0' are phases in the 
hght Majorana mass matrix, one finds that the mixing angles can be 



vr 
4' 



0, 



vr 
4 



(36) 



However, this mixing is not stable under perturbations of the degenerate texture (|3^ ) 
After including renormalisation group effects, the associated mixing angles become 



-0.327, 02 ~ 0.415, 



-0.884 



(37) 



which is inconsistent with a degenerate mass scale much above 1 eV. However, ways around 
this difficulty have been proposed in the context of non-Abelian fiavour symmetries, as 
discussed later. 



4. Neutrino Masses from Flavour Symmetries 
IV- A. Examples based on Abelian Groups 

In this section we give a brief description of models based on extra Abelian symmetries 
which lead to a consistent charged fermion mass spectrum and give predictions for the 
neutrinos and fiavour violating processes. The fact that the fermion mass matrices exhibit 
a hierarchical structure suggests that they are generated by an underlying family symmetry, 
of which the simplest examples are based on Abelian groups. To review how the various 
terms in the mass matrices arise in such a model, we first denote the charges of the Standard 
Model fields under the symmetry as in Table 1. The Higgs charges are chosen so that the 
terms fsf^H (where / denotes a fermion and H denotes Hi or H2) have zero charge. 
Thus, when the U{1) symmetry is unbroken only the (3,3) elements of the associated mass 
matrices will be non-zero. When the U{1) symmetry is spontaneously broken via standard 
model singlet fields, 9, 9, with opposite U{1) charge and equal vevs (vacuum expectation 
values), the remaining entries are generated in a hierarchical manner. The suppression 
factor for each entry depends on the family charge: the higher the net U{1) charge of a 
term fifjH, the higher the power n in a non-renormalizable term fifjH (^-^^ that has 
zero charge. Here M is a mass scale associated with the mechanism that generates the 
non-renormalizable terms. A common approach communicates symmetry breaking via 
an extension of the 'see-saw' mechanism, mixing light and heavy states, known as the 



Froggatt-Nielsen mechanism pl|] . 



•^Renormalisation-group effects on schemes with neutrino degeneracy have also been discussed in Eol 
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Qi ul di Li ei z/f 


f/(l) 


Pi T'i 



Table 1: Notation for the U{1) charges of the Standard Model fields, where i stands for a 
generation index. 



We discuss the simplest possible scheme, with symmetric mass matrices |^ 0. This 
leads to three viable cases with charges |]26| 



A) : bi = Ci = di 



B) : bi = Ci = di 

C) : bi = Ci = di 
leading to three possible charged-lepton matrices : 

f7 z3 



7 1 



5.io 

2'2' 
(3,0,0) 



(3^ 



Mf oc 




Mf 



oc 



oc 



-5/2 



e 

f1/2 



£5/2 

1 



/ ^ ^3 ^3 



e3 



(39) 



The first two matrices lead to natural lepton hierarchies for e ^ 0.2 and imply large but 
non-maximal lepton mixing. On the other hand, the third matrix leads to maximal mixing 
in the (2-3) sector, but requires an accurate cancellation in order to get the correct ratio 

The neutrino Dirac mass is specified to have the same form as the charged leptons, but 
with a different expansion parameter. Indeed, since neutrinos and up-type quarks (charged 
leptons and down-type quarks) couple to the same Higgs, they should have the same expan- 
sion parameter e(e), where the spread between the up- and down-quark hierarchies requires 



e ^ ^ 0.05 [E2|. Then, 




e 

,5/2 



e 

.1/2 



£5/2 

el/2 

1 



(40) 



for the first two choices of charges in (^) respectively. 

The mass structure of the light neutrinos is more complicated, due to the heavy Majo- 
rana masses of the right-handed components. However, as discussed above, the neutrino 

*The lepton sector in this case is identical to that of SU{3)c x SU{3)l x SU{3)r |2^, which, however, 
predicts asymmetric quark mass matrices. 
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mixing in the case of large neutrino hierarchies and non-zero determinant solutions is de- 
termined entirely by the Dirac mass matrices (which feel the left-handed charges). The 
heavy Majorana mass sector only affects the neutrino eigenvalues. For instance, for the set 
B of solutions, with hi = Ci = di = (^|, |, 0^, we find that [26, J3 



-e2 1 

g5/2 _^l/2 



_e5/2 
el/2 



e 
1 



(41) 



whilst for the set A of solutions with 6,- = c,- = di 



.1 I 

2 ' 2 ' 



we find that 



/ 


1 


e2 


_e7/2 






1 


el/2 


V 


^7/2 


-6-1/2 


1 



/ 1 



— e 

z7 



1 

— e 



— e 
e 
1 



7 



(42) 



We see that both of these solutions are close to the limit of 'matched' mixing discussed in 
the previous section, since the mixing in the (1-2) and (1-3) sectors is much smaller than 
in the (2-3) sector, even though the latter may be less than maximal. 

The solution C is also close to the limit of 'matched' mixing, but with maximal (2-3) 
mixing. The lepton mixing matrix in this case is 



_ 1 

2 



1 
2 

1 
2 



1 

V2 I 



(43) 



and Vyj^ has a similar form, with the small entry — > e*, in analogy with (^ ^). The 
diagonalization of the lepton mass matrix leads in this case to eigenvalues e^, — e^, 2, and 
some fine-tuning would be needed to obtain correct low-energy masses. In this case, 
might differ from (^3]). 



IV-B. Examples based on non-Abelian Groups 

We do not discuss these type of models in detail, but remark that they lead naturally to 
models with degenerate neutrinos. Indeed, when the three lepton doublets form a real 
irreducible representation of some non-Abelian flavour group, as for instance a triplet of 



5*0(3), which has been extensively studied during the last year ||29[, one expects exact 
neutrino mass degeneracy at zeroth order. The same result can be achieved by discrete 
non-Abelian symmetries |Q and is to be contrasted to the predictions of abelian flavour 



symmetries, which naturally lead to large hierarchies between the various mass entries, as 
reviewed above. However, the non-Abelian symmetry is typically broken by terms related 
to lepton masses. Once the flavour symmetry is broken, the exact neutrino mass degeneracy 
is lifted by small terms and one may be able to reproduce the Super-Kamiokande data. 

Generic neutrino mass textures that originate from such flavour symmetries are of the 



'mismatched-mixing' type, as in the example (|34D mentioned above. It was pointed out p9 
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that such textures are vulnerable to radiative corrections, which may lead to unacceptable 
patterns of masses and mixing angles. The requirement that such a disaster be avoided 
imposes severe constraints on the mixing angles and requires the mixing should be close to 
bi-maximal [^. Thus, one would expect a large mixing angle in the /i — e flavour sector, 
which would tend to generate relatively large large rates for — >• cy, fi 3e and fi ^ e 
conversion, as well as for r ^ fi'y decay. 



IV-C. Flipped-5't/(5) Model derived from String 

We now outline how the above analysis may be extended to a typical grand-unified model 
derived from string. In such a framework, the following features appear generic: (i) non- 
Abelian symmetries are disfavoured, (ii) the U{1) symmetries and charges are specified 
in any given string model, and one generically expects a product of anomalous and non- 
anomalous Abelian groups, rather than a single U{1), (iii) there are many singlet fields 
involved in the mass generation, not just a single pair 6, 9. Once a string model is chosen, 
e.g., by specifying fermionic boundary conditions on the world sheet, then automatically 
the gauge properties of the model and the quantum numbers of all fields, including those 
which may acquire non-zero vevs and fill in the fermion mass matrices, are specified. The 
field vevs that determine the magnitudes of the various entries are constrained by the 
anomaly-cancellation conditions and the fiat directions of the effective potential in the 
theory. Finally, we recall that (iv) additional string symmetries (expressed through calcu- 
lational selection rules |^6|) further constrain the possible forms of the mass matrices, since 
they forbid most of the Yukawa couplings that are allowed by the rest of the symmetries 
of the model. 

An example of the above class of models is provided by the flipped SU (5) model derived 
from string specializing to the pattern of vevs and mass matrices discussed in P7| , 



Looking at the field assignments in group representations, one sees that: (i) since the 
charge conjugates of the right-handed neutrinos have the same charges as the down quarks, 
the Majorana mass matrix will be constrained, and (ii) due again to the common charge 
assignments, the Dirac neutrino mass matrix is the transpose of the up-quark mass matrix. 
The quark and charged-lepton mass matrices have been presented in , where the possible 
fiat directions of the theory were also reconsidered. Since the analysis of the surviving 
couplings after all symmetries and string selection rules are taken into account is quite 
involved, we refer to the original references for details, while here we just give an illustration 
of the predictions for neutrino masses. 

Within this model, we found that the charged-lepton mixing matrix is given by 

Vi= \ -X 1 - I (44) 


where x is the vev of a combination of hidden-sector fields that transform as sextets under 
5*0(6), that needs to be OiV) for realistic quark mass matrices [^. The Dirac neutrino 



mass matrix, m^, is expressed in terms of three expansion parameters f,x,y related to 
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field vevs scaled by the string mass: 



m 



D 



/ x/ 1 
/ X 
V y 



(45) 



where y stands for the SU (5) decuplets that break the gauge group down to the Standard 
Model, with scaled vev ^ Mqut/Ms, where Mg is the string scale. In weakly-coupled string 
constructions, this ratio is suppressed, but the GUT and the string scales can coincide in 
the strong-coupling limit of M theory, in which case y ~ 1. Finally, / stands for a singlet 
field, the value of whose vev ^ 0.04 is fixed in order of magnitude by the quark mass 
hierarchies. 

Since now the mixing in the (2-3) sector of the Dirac mass matrix is zero, we see that in 
order to generate a large hierarchy among the various neutrinos, we necessarily require zero- 
determinant solutions. In where the expectations for neutrino masses were studied, 
we ended up with two possible forms for M^,^, depending on the vevs of the singlet fields. 
These were 



where in the second example a notional factor of A ~ (9(1) has been included so as to avoid 
sub-determinant cancellations, which are not expected to arise once coefficients of order 
unity are properly taken into account. 

As can be seen from the lower 2x2 matrices in (^61) , consistency with the neutrino 
data imphes ?/ ~ 1, as could occur in the strong-coupling limit of M theory, and t f . 
The parameter a represents a higher-order non-renormalizable contribution , where 

$ is an effective singlet, that is expected to appear at some power n > 7. Its actual value 
is irrelevant, provided that it not too tiny, in which case m^-^ might be increased to an 
unacceptable value. 

The situation described just above corresponds to the 'mismatched' scenario introduced 
above and exemplified previously in the context of simpler U{1) models. It is straightfor- 
ward now to determine the effective light neutrino mass matrix for the above cases, and see 
that large mixing is implied for the — z/^ sector, as required by the atmospheric neutrino 
data. We do not go into further details, since they would depend on more specific aspects 
of the model, which contains several poorly-constrained expansion parameters, not all of 
which are necessarily very small. However, we can infer some qualitative properties of the 
predictions of this fiipped SU (5) model for fiavour- violating decays. From the form of the 
charged- lepton mass matrix (0), we would expect a rather large /i — 67 branching ratio. 
On the other hand, the r lepton remains completely decoupled from the lighter families in 
the approximation considered here 0. This is in contrast with many of the U{1) models, 
with the result this fiipped SU{5) model would predict a relatively small branching ratio 
for r —>■ fi^y. 

^This is, of course, connected to the fact that the large v^ — Vt mixing needed to interpret the atmospheric 
neutrino data comes entirely from the effective light Majorana matrix. 




M^^ oc Xxy'^ 
V fy 




(46) 
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5. Predictions for Rare Processes 



We now discuss quantitatively implications of the large neutrino mixing needed to inter- 
pret the neutrino data for processes that violate the conservation of charged lepton flavours. 
In particular, as we shall see, the likelihood that the atmospheric neutrino problem is solved 
by z/^ — v-r mixing suggests the likely appearance of the t —>■ fij transition at a rate that 



may be accessible fl^ . This observation supplements the common belief that fi ^ ej and 
related processes may offer good prospects for observing charged-lepton-fiavour violation. 

In this Section we present results for these and related flavour-changing decays for 
generic examples of the U{1) textures discussed in the previous Sections. We start with 
the radiative decays /i — 67 and t ^ fi'j 0. For the reasons discussed earlier, we do not 
explore further the corresponding decays where the photon is replaced by an e"'"e~ pair. 
However, we do present later some numerical results for fi — e conversion on Ti, whose 
rate is not directly related to that for n —>■ ej, and which may present some experimental 
advantages |jlO| . 



In addition to the flavour-mixing effects, the rates for flavour non-conserving decays 
are also sensitive to other physical quantities. In particular, in the supersymmetric GUT 
context explored here, they depend on the masses of the sparticles that mediate the flavour 
non- conserving processes. We parametrize their masses in terms of the universal GUT-scale 
parameters mo and mi/2, and use the renormalizat ion-group equations of the MSSM to 
calculate the low-energy sparticle masses, taking into account low-energy threshold effects. 
Other relevant free parameters of the MSSM (discussed in detail in 0) are the trilinear 
coupling A, the sign of the Higgs mixing parameter /i, and the value of tan (3. Here we fix the 
value of Aq = —mi/2, assume the sign of the parameter to be either positive or negative, 
and restrict our analysis to low and intermediate values of tan/5 < 10. Before presenting 
our results, we first consider in more detail the relevance of some of the parameters entering 
in the calculation. 

• The case of non- universal soft masses at the GUT scale was analysed in ^ . However, 
in such a case, the predictions for fi —>■ ej exceed the current limits for most of the 
parameter space, this is why, in the present work, we restrict ourselves to the case of 
universal soft masses at the GUT scale, assuming the existence of some mechanism that 
assures universality and the absence of mixing in the Kahler potential. 

• No constraints on the mo,mi/2 values were discussed in [Q. Here, however, we choose 
their initial values so as to respect the cosmological relic-density constraints, as discussed 
in 



m 



• To calculate within a given fermion mass texture, we need to estimate the mixing 
both of the neutrinos and of the charged leptons. In the case of U{1) models, all the mass 
matrices are given as power expansions in the parameters e, e, and it is not possible to fix 
uniquely the C(l) coefficients Cij needed to obtain masses and mixings. There is a similar 
ambiguity in flipped SU (5) models. 

For our numerical study here, we adopt for the three U{1) models discussed in Section 4 
an approach similar to that of [S], commenting later on the expectations for other classes of 



models. In f^, a set of 0{1) coefficients was used which led to the correct mass spectrum, 
with the mixing in the charged-lepton mass matrix emerging as a prediction. Nevertheless, 
the choice of Cij is not unique, since many sets can give the same mass eigenvalues but 
different mixing. In order to test the sensitivity of our results to this arbitrariness, here we 
use three representative sample sets of coefficients Cij = Cji for each of the U{1) textures 
presented in Section 4, namely 

A) ci2 = .87, C22 = .59; C12 = .89, C22 = 1.44; C12 = .61, C23 = .79 

B) C12 = 1.68, C22 = .56; Cu = 1.70, C22 = 1.45; Cu = 1-41, C23 = .79 

C) C12 = 2.59, C22 = 1.27; cu = 2.34, C22 = .79 cu = 2.75, C23 = 1-13 

which we denote hereafter by ^1,2,35-81,2,3 and Ci,2,3, respectively. 

We plot in Figs. 2, 3 and 4 the branching ratios for /i — >• 67 and r — > /i7 in the textures 
A, B and C, respectively, assuming in each case the representative value mi/2 = 250 GeV, 
tan/3 = 3 and allowing both possible signs of /x. In each of Figs. 2, 3 and 4 the spread in 
the numerical results reflects the uncertainties associated with the numerical coefficients 
Cij. This spread immediately warns us that the results shown in Figs. 2, 3 and 4 should 
be understood only as order-of-magnitude estimates. We note that the predictions for 
the /i — > 67 decay branching ratio in textures A and B are particularly sensitive to these 
coefficients, whilst r —> /i7 reaction is generally less sensitive, as also is fi —>■ ej in texture 
C. In all cases, the /i > choices (solid curves) lead to an enhancement of the branching 
ratios compared to the choice /i < 0, at least for larger values of mo. We note also that, 
for a small range of mo values that depends on the texture and the specific choice of the 
coefficients, there is potentially a strong cancellation of the mixing effects, leading to a 
considerable suppression of the branching ratio. 

Among the three cases studied, texture B generally gives the highest predictions for 
BR{^ — ^ 67). This is, of course, a simple consequence of the fact that the (1,3) mixing 
matrix element is larger in case B than in case A. Texture C, on the other hand, predicts 
BR{fi — * 67) < 10~^^ for all Cij choices. As for BR{t fi'j), we note that it is in general 
enhanced for small mo values. Also, texture A may yield BR{t fi'y) as high as 10"''. 

We present predictions for the larger value tan/3 = 10 in Fig. 5, for both the fi —>■ ej 
and r — >■ /i7 branching ratios, using the three sets of coefficients Cij for texture A. Again, 
the plots are obtained using the representative value mi/2 = 250 GeV and the two possible 
signs of yU-parameter. The results are qualitatively similar to those for the tan (3 = 3 case 
in Fig. 2, though the fi ^ e'y branching ratio now tends to lie above or close to the present 
experimental limit for two choices of the coefficients Cij. Only in the third choice of Cij is 
BR{fi 67 well below the present experimental limits. We also observe in the right-hand 
plot of Fig. 5 that the branching ratios for r — > /i7 are enhanced compared to the low tan /3 
case. 

We show in Figs. 6, ^predictions of the branching ratios for the radiative decays fi —>■ 
67, r fi'y in the (m.o,mi/2) plane, specializing to texture A, using tan/3 = 3 and 10 and 
assuming /i < 0. In general, the branching ratios tend to decrease as mi/2 increases. If 
tan /5 = 3 as shown in Fig. H, case Ai predicts values of BR{n 67) compatible with 



-17- 



10-" 



10 



A 

I 

£'io " 

CO 



10 



10-'= 





' 1 1 1 

Experimental bound 




~ — 

A3 




A3 











100 



200 



300 



400 



500 



10° 
10-= 
10-' 
^10-= 

=1 

^10-= 

10-" 
10-" 
10-" 




100 



200 



300 



400 



500 



m„ 



m. 



Figure 2: Predictions for BR{^ 67) and BR{t — > /I'j) for texture A of Section 4, O'S- 
suming the values mi/2 = 250 GeV, tan/? = 3 and Aq = —mi/2- The solid lines correspond 
to positive fi, and the dashed ones to fi < 0. The results are for the three specific choice of 
the undetermined numerical coefficients Cij shown in the text. We see that, for fixed uiq, 
the II ^ curves are more sensitive to the Cij than are those for t ^ iij. 




Figure 3: As in Fig. 1, hut for texture B of Section 4- 
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Figure 4: As in Fig. 1, hut for texture C of Section 4- 




Figure 5: As in Fig. 1, hut for tan (5 — 10. 
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Figure 6: Contour plots in the {'mi/2,mQ) plane for the decays fi ^ cj (left) and t —>■ fi'j 
(right), assuming tan (3 = 3 and /i < 0, for the cases ^41,^3. We see that the rates for 
both decays are encouraging throughout the dark-shaded region preferred by astrophysics 
and cosmology 



the experimental bound in most of the cosmologically preferred region. In contrast, if 
tan /5 = 10 as shown in Fig. acceptable BR{fi — >• ej) rates are found only for large values 
of mo > 400 GeV). In this latter case, the choice of coefficients is more favoured. We do 
not display similar plots for the other two textures, but note that texture B is more sensitive 
to the choice of coefficients Cij, and that for certain choices of them it predicts values close 
to the experimental bounds for a large portion of the plane mo, mi/2 < 500 GeV. 

The light-shaded areas in Fig. ^ correspond to the regions of the {mi/2,mQ) plane 
that are excluded by LEP searches for charginos and by the requirement that the lightest 
supersymmetric particle not be charged [Q. The dark-shaded areas in Fig. ^ are those 
where the cosmological relic density is in the range preferred by astrophysics. We see that 
both the decay modes n —>■ cy and t ^ fi'y could well be measurable throughout this 
astrophysical region. In particular, we find BR{t — > fij) > 10~^ in most of this region, 
and for mi/2 < 220 GeV and mo < 110 GeV there are regions where BR{t — * fi'-f) > 10"^. 
These observations also apply to texture B, and to the choice tan/3 = 10 (not shown). The 
predictions of texture C may also reach above 10~^, but reach above 10~® only in a small 
portion of the cosmologically-favoured region. Texture B leads to the highest predictions 
for BR{fi —>■ cy), typically above 10~^^ for most of the cosmologically-favoured region of the 
(mo, mi/2) plane, and even reaching values above 10~^^ in a small portion of the parameter 
space. The sensitivity of the branching ratios to different choices of coefficients Cij can 
be seen by comparing the solid and dashed lines in Fig. ^. For tan/5 = 3, the obtained 
values for the branching ratios under consideration are smaller. A considerable sector of 
the cosmologically-favoured region, leads, however, to values for the branching ratios of the 
same order as the ones discussed above. 

We show in Fig. 8 the correlation between the n —>■ cj and r —>■ fi'j branching ratios 
for two different choices of coefficients Cij for each of the textures A, B and C presented 




Figure 7: Contour plots in the (^7^1/2, "^o) plane for the decays fi —>■ cj (left) and r —>■ 
/X7 (right), assuming tan/5 = 10 and /i < 0, and using the cases Ai (solid lines), A3 
(dashed). The rates for both decays are enhanced relative to the small tan/3 case. As a 
result, case A3 is more favourable now in the dark-shaded region preferred by astrophysics 
and cosmology 



in Section 4. In each case, the branching ratios have been calculated for a sampling of 
(mo, mi/2) pairs in the cosmologically-favoured region of |^ for tan f3 = 10: similar results 
hold for tan (3 = 3. We see clearly how the correlations between the two branching ratios 
vary with the choices of the Cij coefficients, because of their influences on flavor mixing. In 
the two cases Ai and ^3 (texture A2 is similar to Ai) presented in Fig. 8, the values of 
the predicted branching ratios are characteristic for each Cij set, with the case ^3 generally 
predicting smaller ratios for BR{fi e'j)/BR{T yU7). In the case of texture C, the 
dependence of the results on the choice of Cij coefficients is rather different: this texture 
tends to predict a relatively large ratio BR{fi cy) / BR{t — > /i7) . 

Finally, we show results for /i — e conversion, using both penguin and box diagrams. We 
gave in Section 2 an order-of-magnitude estimate of the branching ratio for this reaction, 
namely BR{^Ti eTi) ^ 5.6 x 10^^ BR{fi — 67). However, as we commented there, the 
two processes exhibit different functional dependences on the A'^^^ functions (indeed, as 
is well known, /i — >■ 67 does not depend at all on the parameters r). Further, in the 
case of /i — e conversion there are additional contributions from box graphs which further 
complicate the ratio BR{fiTi — >■ eTi) / BR{n — > 67) . In view of the possible improvement of 
the experimental sensitivity to fi — e conversion, we present plots similar to given previously 
for /i — > 67. 

It is instructive to compare the fi — e conversion rate with the corresponding predictions 
for BR{fi 67), to explore the effects of penguin and box diagrams. In Fig. 9 we plot 
the ratio of the /x — e conversion and /i — 67 rates versus the scalar mass parameter mo. 
We see immediately that the ratio is not constant, though its general order of magnitude 
is that estimated in (pOI). The dependence of the penguin contribution on mg is shown 
separately from the combined effects of penguin and box diagrams, for the two signs of /x. 
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Figure 8: Scatter plot of model predictions for BR{n — > ej) versus BR{t — >• H'^), assuming 
tanj3 — 10 and // < 0. The circles correspond to texture Ai, the + to texture A2, the stars 
to texture Bi, the x to texture B2 the squares to texture Ci, the diamonds to texture Ci- 
Note the characteristic correlations in the different models. 

In the case of > 0, BR{ii — > 67) is relatively enhanced for large mo values, whilst the 

opposite is true for fi < 0. 

We show in Fig. 10 the dependence of the //. — > e conversion rate on mo for the two 
textures A, B and for the values tan (3 = 3 and mi/2 = 250 GeV. As before, we choose three 
sets of numerical coefficients Cij for each texture, and exhibit the results for both signs of 
Interestingly, cases Ai, A2 give a rate close to the present experimental limit for most of the 
mo region explored when fi > 0. The corresponding predictions for < are considerably 
lower for large mo values, but converge with those of the fi > case for small mo- We note 
that texture B exhibits greater sensitivity to the coefficients Cij than docs texture A. 

Fig. 11 displays contours of the rate for — > e conversion in the (mo,mi/2) plane for 
the two scenarios Ai and ^43. We see that the former predicts a rather larger rate, which 
offers good prospects for observation throughout the region preferred by cosmology. On 
the other hand, scenario A3 predicts a rather lower rate for // — > e conversion. 

We comment at this stage on some of the ambiguities in the results shown above, 
within the general U{1) framework studied. We have already explored to some extent the 
ambiguity associated with 0{1) coefficients in the Dirac mass matrices for the fermions. The 
ambiguity induced by simple sign changes can be particularly acute, as we illustrate with 
one simple exercise. For example, since the light Majorana mixing matrix is Vmns — V^Vuj 
any modification that changes the relative signs of off-diagonal entries in and could 
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Figure 9: The ratio of the /i — > e conversion rate to BR{fi — > 67) is plotted versus tuq, for 
H < 0. We plot the penguin contribution separately from the sum of the penguin and box 
diagrams, /or tan /3 = 3 (left) and tan/3 = 10 (right). 

cause large changes in the mixing angles of Vmns, as one changes from destructive to 
constructive interference, or vice versa, with intermediate possibilities corresponding to 
various phase possibilities that are not specified by the U{1) symmetry. 

As an exercise, using the numerical values of the coefficients in the case Ai discussed 
above, we invert the signs of all the off-diagonal entries in m^, and repeat consistently all 
the subsequent steps in the calculations. The implications for /i — > 67 and r — > /i7 are 
shown in Fig. |T2|, where we see that BR{fi cy) may be increased by about two orders of 
magnitude 0, whereas that for BR{t fi'y) is increased by less than an order of magnitude. 
This 'inverted' U{1) model actually has much larger — Uj. mixing than the 'uninverted' 
version of Ai, or indeed the other textures studied previously, agreeing better with the 
Super-Kamiokande data. We interpret the difference between Ai and 'inverted' Ai, on 
the one side, and between Ai and ^3, on the other side, as indicative of the numerical 
ambiguity within any particular texture [[ We see that the predictions for — > 67 are very 
variable, bracketing the present experimental upper limit within a broad range, whereas 
the predictions for t ^ fi'y are more closely bunched, and likely to be within reach of 
experiment. 

The generality of these features is left for exploration in the future. There are also 
other ambiguities, for example in the choice of the heavy Majorana mass matrix 
whose detailed study we also leave for future work. 

Before concluding, we make a few more remarks about non-Abelian models, and about 

®We find a similar enhancement for the /i — > e conversion rate. 

^The drastic changes in the 'inverted' Ai case may not be possible in GUT models where the neutrino 
Dirac matrix is related closely to the up-quark mass matrix. In such models, large — i^r may also be 
arranged by a suitable choice of the heavy Majorana sector. 
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Figure 10: The muon conversion rate ji ^ e as a function of for the textures A and B, 
assuming tan/3 = 3. The solid lines correspond to /i > and the dashed to /i < 0. 

the flipped SU{5) model. As commented at the end of subsection IV-B, generic non- 
Abelian models would fall into our 'mismatched' category, and we would expect them to 
have relatively large rates for both — e and r — > // transitions, as a result of their 
preference for near-bi-maximal mixing. Therefore, we expect the prospects for charged- 
lepton-flavour violation in these models to be at least as favourable as in the Abelian 
models studied here. In the case of flipped SU{5), if we use naively the matrices displayed 
in Section 5, we find that the // — > 67 process is rather enhanced, and it exceeds the 
present experimental bounds in a considerable region of the parameter space, at least for 
some generic choices of undetermined numerical coefficients. However, since the plethora 
of poorly-constrained expansion parameters introduces ambiguities, as discussed earlier, a 
complete exploration of this model is beyond the scope of this paper, and it may be that the 
model can survive for suitable values of these coefficients. On the other hand, the r ^ /X7 
reaction is highly suppressed, because in this model all the mixing needed to interpret the 
atmospheric neutrino data comes from the neutrino sector, and there is no mixing in the 
T — 11 charged-lepton sector. 



6. Conclusions 

Although family symmetries provide many interesting insights into the hierarchy of 
fermion masses, there is no unique framework that fits the available information on charged 
fermions and neutrinos. As a result, there is considerable ambiguity, even within the 
subclass of Abelian flavour models, in their predictions for charged-lepton-flavour violation. 
We have explored some of the range of possibilities in this paper. Many sets of undetermined 
0{1) coefficients in Abelian flavour symmetry models of the charged lepton mass matrix 
can fit well the mass spectrum and the neutrino data, but vary in their predictions for 
the branching ratios of rare processes. As examples, three different sets of coefficients 
were used to fit the charged lepton and quark mass matrices in each of three Abelian 
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Figure 11: Contour plots in the (mi/2, '^o) plane for /i — > e conversion, assuming tan /3 = 3 
and /i < and using the cases Ai,A3. We see that the conversion rate is encouraging 
throughout the dark-shaded region preferred by astrophysics and cosmology in the sce- 
nario Ai. 

texture models. Our studies show that these vary in their predictions for flavour-changing 
branching ratios by up to two orders of magnitude. The good news, however, is that many 
of these models seem to be accessible to a new round of experiments searching for /i — e 
and r — /i transitions. In particular, we would like to re-emphasize the interest of exploring 
the branching ratio for t ^ fi'y down to the 10~^ level or below, as may be possible at the 
LHC. 

In view of their accessibility, and precisely because of their model-dependence, such rare 
decays may become a powerful tool for distinguishing between different neutrino textures. 
As we have seen, textures based on Abelian flavour symmetries tend to predict relatively 
small /i — e flavour mixing, thus leading to rates for /i — > 67, /x — 3e and /i — e conversion 
in nuclei that are generically below the experimental bounds, though close enough to offer 
interesting physics opportunities for experiments with present and future intense /i sources. 
On the other hand, schemes based on non- Abelian flavour symmetries would tend to predict 
large mixing in the (1-2) lepton sector. In this case, larger rates are likely to be found for 
the above processes, and in certain textures part of the supersymmetric parameter space 
may already be excluded. It is likely also that string-derived flipped SU{5) schemes based 
on would have large off-diagonal entries in the (1-2) lepton sector, leading to larger rates 
for /i —>■ cj, fi 3e and /i — > e conversion in nuclei than the simplest models with Abelian 
flavour symmetries, though the rates for r — ^^7 would be relatively low. 

We conclude by encouraging the community to pursue actively new generations of ex- 
periments to probe charged-lepton-flavour violation, in both fi and r decays. Such efforts 




Figure 12: Predictions of BR{fi 67) and BR{t — ^ fi'y) in the 'inverted' version of Ai, 
obtained by inverting the signs of all the off-diagonal terms in , for fi < 0. 

would complement nicely the physics being revealed by neutrino oscillations, and could 
provide precious insight into the flavour problem. 
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Note Added 

While this paper was in the flnal stages of preparation, we received the paper by Feng, Nir 
and Shadmi in 0, which makes points similar to ours, in a complimentary way. 
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